Pion form factors 
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We consider the electromagnetic and transition pion form factors. Using dispersion relations 
we simultaneously describe both the hadronic, time-like region and the asymptotic region of large 
energy-momentum transfer. For the latter we propose a novel mechanism of Regge fermion exchange. 
We find that hadronic contributions dominate form factors at all currently available energies. 



o 

(N 

(N 
I 



> 
oo 
in 
in 
in 

o 



> 
X 

S-H 



Photons interact with quarks, the charged con- 
stituents of hadrons and the resulting electromag- 
netic form factors probe the quark energy-momentum 
distribution in hadrons. In this letter we ex- 
amine the charged pion electromagnetic form fac- 
tor F2tt{s), which is defined by the matrix element 
(7r+(p')7r-(p)|J^|0) = e(p' - p)fj,F2T,{s), and the tran- 
sition from factor between the neutral pion and a real 
photon, Ft^^{s) determined by (7r°(p')7(A,p)| J^JO) = 
ie^/47r^/^e^,.a/3e*'^(A)p'"p'^i^7r7(s)- Above, is the elec- 
tromagnetic current, s = {p' +pY is the four-momentum 
transfer squared and /tt — 92.4 MeV is the pion decay 
constant. Current conservation implies i^27r(0) — 1 and, 
in the chiral limit, axial anomaly determination of the 
7r° — > 27 decay leads to the expectation Ft^^{Q) w 1. 
Because at short distances quark/gluon interactions are 
asymptotically free, it has been postulated that at high 
energy or momentum transfer, \s\ 3> /i^, both form fac- 
tors measure hard scattering of the photon with a small 
number of the QCD constituents [TH3]. One would then 
expect ~ 0(1 GeV^), which is the typical hadronic 
scale, however, given the current status of the data it 
seems that ^2 could be as large as ©(lO-lOOGeV^) [1[5]. 
This implies that an alternative description of the un- 
derlying dynamics might be in order and the subject of 
applicability of pQCD to exclusive reactions has in fact 
a long history [5]. Perturbative QCD (pQCD) analy- 
sis of the form factor asymptotics assumes specific prop- 
erties of certain non-perturbative quantities, i.e. the 
parton momentum distribution amplitudes in the low- 
momentum," wee" region. If these had different behavior 
from what is assumed in the pQCD analysis the argu- 
ments leading to dominance of leading twist perturba- 
tive scattering would break down [7] . Such pion distribu- 
tion amplitudes were considered recently in [H [S] , how- 
ever, the authors used perturbative evolution to soften 
the wee region and use pQCD formulae. The available 
data on the pion electromagnetic form factor ranges up to 
|s| < 10 GeV 10 and is approximately a factor of three 
above the asymptotic prediction ^llj . Even more spectac- 
ular discrepancy is observed in the transition form factor 
recently measured by BaBar [12] . For momentum trans- 
fers as large as — s « 40GeV^ the measurement disagrees 
with the asymptotic prediction not only in normalization 
but also in the overall s-dependence. While pQCD pre- 
dicts sFTr-Y(s) 2/^ as \s\ —7- 00 [3], the data suggest 



that the magnitude of ^sFT,~f{s) grows with \s\. 

Crossing symmetry implies that form factors in the 
space-like (s < 0) and time-like (s > 0) region are 
boundary values of an analytical function defined in the 
complex-s plane with a unitarity cut running over the 
positive s-axis and starting at the two pion production 
threshold branch point Sth = 4mJ. In the time-like 
region the electromagnetic (transition) form factor de- 
scribes the amplitude for production of a spin-1, 7r"'"7r~ 
pair (7r''7) in the external electromagnetic field of the vir- 
tual photon. In the space-like region the form factors are 
usually interpreted in terms of parton three-momentum 
distribution in a hadron (and/or photon). Analyticity de- 
mands these apparently distinct physical pictures to be 
smoothly connected. The dominant feature of the spin- 
1, 7r+7r~ state is the isovector p(770) resonance, which 
also dominates the electromagnetic form factor. There is 
no time-like data available for the transition form factor, 
however, also in this case one expects to see the p, and 
the isoscalar, w(782) resonance. The analytical continua- 
tion to the space-region implies that for —s < IGeV^, i.e. 
in the hadronic range, the quark wave function is dual to 
the vector-meson exchange in the crossed channel. 

In the following, we relate the space-like and time-like 
regions through dispersion relations, and focus on the dy- 
namics in the asymptotic region, s — ?> +00. In view of the 
BaBar "anomaly" and the apparent failure of the pQCD 
description, we propose a novel description for the dom- 
inant mechanism that drives the asymptotic behavior of 
the form factors. 

The discontinuity of Fj^-y^s) across the unitary cut is 
given by 

X 

. (1) 

and the sum runs over all possible intermediate states 
X 7^ 27r, Stt. Here, tx^rr-y (Fx) represent the amplitudes 
for X — > 7r''7 (7* — >• X), respectively and px is a prod- 
uct of the phase space and kinematical factors {i.e. for 
the 27r intermediate state P2-k{s) — s{l — sth / s)'^^^ /96t:). 
Provided ImF^^ vanishes at s — )■ 00, its real part can be 
reconstructed for any s from the unsubtracted dispersion 
relation 



F^j (s) 



ds 



,ImF^-f{s') 



(2) 



2 



The two lowest mass intermediate states, X = 27r, Stt 
that are dominated by the p(770) and w(782) resonances, 
respectively, are expected to saturate the cut in the 
hadronic range sth < s < IGeV^. The a;(782) in the 
isoscalar Stt channel is a narrow resonance with width to 
mass ratio, T^^/mi^ ^ 10^^ and its contribution to 
can be well approximated by a Breit-Wigner distribution. 



with c^y"^ 



TO' 



im„r„(s) 



(3) 



0.493 obtained from 



uj ^ TTj and oj — )■ e+e~ decay widths yielding gi^-rr-i — 
1.81 and g^^ = 17.1, respectively. The contribution from 
the 2tt intermediate state is dominated by the p(770) 
resonance, which determines both the t2;T,7r7 scattering 
amplitude and the pion electromagnetic form factor, i^2ir 
for s < IGeV^, and vector meson dominance (VMD), 

yields c^j'^ = 47r^/^5p7r7 /wp^p = 0.613 (with p — 717 
and p — >■ e"'"e decay widths leading to gpT^^y = 0.647 and 
gp = 4.96). At s = the sum of the two resonance con- 
tributions to Fj^j agree with the anomaly driven normal- 
ization to within 10-15% and the isovector contribution 
can be further improved using a unitary parametrization 
of [13l[T4j, which for i^2ir and i27r,7r7 in Eq. ([T]) yields, 

F^As) ^ PisMs), i2.,.^(s) = C-i(s)l](s), (4) 

where ri(s), the Omnes-Muskelishvilli function [TS] com- 
puted from the phase of the vector-isovector elastic rnr 
scattering amplitude and satisfying the VMD relation, 
n(s ~ TO^) ~ ■m?p/{ml - s - impVp{s)). The poly- 
nomials P(s) and C{s) {P{s) = 1 + O.Is/to^, C{s) = 
fl[l + 1.21 si ml + IMs^/jTv-p - O.SOsVm^]) are deter- 
mined from fits to the electromagnetic form factor and 
the solution to the dispersion relation for the t2Tr,7r7 am- 
plitude, respectively. At higher energies, s > IGeV"^ the 
KK inelastic channel and other multi-particle interme- 
diate states are expected to contribute. Unfortunately, 
since no time-like data are available (unlike in the case of 
the electromagnetic form factor) one cannot unambigu- 
ously determine these contributions. A possible deter- 
mination of the multi-particle hadronic states could be 
given in terms of quark/gluon intermediate states, much 
like in the derivation of QCD sum rules (c/. Ref. pK^ for 
the case of the pion electromagnetic form factor). 

Since the electromagnetic form factor Fx of a com- 
posite state decreases with energy-momentum transfer, 
asymptotically the r.h.s of Eq. ([T]) is dominated by the 
X — qq, quark-antiquark intermediate state. Its form 
factor is Fqq = 1, (in units of the quark charge) and the 
state contributes to ImFj^^ via the qq — >■ ttj, P-wave 
scattering amplitude, igg,,r7 as illustrated by the last dia- 
gram in Figjl^. The qq contribution shown in Fig|l^ may 
be compared to the one in FigjljD, which represents the 
asymptotic contribution as predicted by pQCD. In the 
latter, the qq — > scattering amplitude, shown to the 



right of the vertical cut line, is given by a free quark prop- 
agator exchanged between the final state pion and pho- 
ton. In the kinematically relevant domain s ^ t ^ b^^ 
with t being the four momentum squared carried by the 
exchanged quark and b « few GeV^^ the typical slope 
of the product of residual coupling of the exchange quark 
(/3^, /3^), the amplitude tqq .^^ is expected to have a Regge 
behavior [T^ 



V-^^(s,0 = /3^(i)/3^(t)s"'(*) « e 



(5) 



The difference between the free, Fig[T]3 and the Regge 




FIG. 1: Hadronic and asymptotic contributions to the tt" 
transition form factor. 

propagator Fig{T^ can originate from the sum of ladder 
gluons in the wee region (c/. Vig^). The quark Regge 
trajectory aq{t) « Oiq{0) -\- a'qt is not known, however, 
phenomenologically it can be related to the leading Regge 
exchange in tttt scattering, i.e. the p or /2 exchange (that 
are nearly degenerate, cf. [H]), ap{t) ^ afrj{t), 



aq{t) - 0.5ap(t) -I- 0.5 w 0.75 + 0.45i. 



(6) 



It is worth noting that the dominance of quark-exchange 
(or, more precisely, quark interchange) mechanism has 
previously been observed in hard scattering processes 
with hadrons TW]. Hard scattering data are furthermore 
compatible with aq{t = -1 GeV^) sa 0.3-0.4 |20j, which 
is consistent with Eq. ^ . Detailed derivation of Eq. ^ 
i.e. relation between quark and meson Regge trajectories 
will be given in the forthcoming paper |21j . After pro- 
jecting onto spin-1 partial wave, the energy dependence 
of the asymptotic, qq contribution to ImF^T^ is therefore 
expected to behave as (modulo terms ~ O(lns)), 



ImF, 



(qq) 



(s) c 



(99)„Q,(0)-3/2 



2ir 



(7) 



The important point is that with aq{0) = 1/2 + e (assum- 
ing a'q{0) > 0, it is consistent with the absence of a physi- 
cal pole for a confined quark) Eq. ([t]) implies asymptotic 
increase of the energy weighted transition form factor. 
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sF^^{s) oc s*^ in agreement with the BaBar measurement. 
Such an increase cannot be accounted for by the exchange 
of the free quark as it is the case for the leading twist 
pQCD model. Combining the uj and the p resonance 
contributions of Eqs. (|3|,(|4| with the asymptotic form 
of Eq. Q and making the simplifying assumption that 
the first two contribute to ImF^^ for s < 1 GeV^ while 
the asymptotic part saturates ImF-^j for s > /i^, we fit 
the the available data using Eq. ^ with the single free 
parameter ci?-f that determines the normalization of the 
asymptotic contribution. The result is shown in Figj2j It 
is worth noting that even at largest values of — s the bare 
qq production gives only about 50% of the form factor 
with the remaining half coming from the resonances. 
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FIG. 2: Our results for 1-F'7r7(s)l in the space-like region for 
p^ = 1 GeV^ (dashed line), 5 GeV^ (solid line), 10 GeV^ 
(dotted line), in comparison with the experimental data from 

[12 [221. 



In the case of the pion electromagnetic form factor, the 
discontinuity across the unitary cut is given by 

X 

(8) 

where in the sum is over intermediate states {X ^ 
2tt, KK) in 7* — > X — > 27r and where we separated 
the two channels X = 27r and X = KK which phe- 
nomenologically are most significant in the hadronic do- 
main |23j . Above the inelastic threshold, s > s^, the uni- 
tarity relation now involves both ImF2-n and ReF2Tr and 
can be solved algebraically. Assuming that the elastic 
amplitude, t27r,27r asymptotically approaches the diffrac- 
tive limit, t27r,27r ^ */2p27ri from Eq. (Isl) one finds 



F2AS) 



2i 



X^2Tr 



tx,27rPxFx 2itqq^2-K OC is 



a,(0)-3/2 



(9) 

Except for the expected energy dependence c/. Eq. ([5]), 
we do not know separately the real and imaginary parts 
of tqq^2-n. Assumlug, as in the case of the transition 



form factor, that the real part of the discontinuity due 
to qq state has the energy dependence given by the reg- 
gized quark exchange, we can compute F2-K using Eq. Q 
and the Cauchy representation (the imaginary part of 
t*q 2TrPx would then be given by the solution of an alge- 
braic equation that follows from Eq. ([8|). This yields, 
F2^{s)=N{s)/D{s) with 



Nis) 
Dis) 



E- 

X^27T 



ds 



~ cxp 



, D{s')Re[t*x,2As')px{s')Fx{s')] 

[l-*^k2.(s'W(sO](s'-s) 



(10) 



The phase is obtained from the elastic amplitude, 
tan(/) = Rei27r.27r/327r/(l — Iiiif27r.27r/027r)- As discusscd ear- 
lier, the dominant feature of the pion electromagnetic 
form factor is the p(770) resonance. Close to the reso- 
nance peak there is also a contribution from the isospin- 
violating a; — ?> 27r decay. Here we do not attempt to re- 
produce the details of the p — uj interference region. The 
next relevant feature is the large variation in magnitude 
of |-F27r| 'At ^y s ^ 1.7 GcV which is also seen in the elastic 
2tt — )■ 27r amplitude and is attributed to the contribution 
from the inelastic p"(1700) resonance decaying to KK. 
We thus approximate the sum over inelastic channels in 
Eq. ( 10 1 by the single KK channel, and above s > p? 



the residual qq continuum with 



R(^i*qq,2-^Px 



- ^(9'3)o",(0)-3/2 



(11) 



For the t 



27r,27r 



and t 



KK,2tt 



amplitudes we use the 
parametrization from |26j . Even though the contribu- 
tion to the dispersive integral from the high energy tails 
of t^^ 9^ and tj^j^ 2-K are suppressed by form factors c/. 



Eq. (|8| we nevertheless extend the parametrization from 
P5j to higher energies by smoothly joining the resonance 
region to the spin-1 projected Regge limit of tttt — >■ tttt 
and KK — >■ tttt amplitudes. We parametrize the isovec- 
tor kaon form factor Fx using Breit-Wigner distribu- 
tions which include the p(770), p'(1400) and p"(1700) 
[24j. Finally we fit the available data on |F2t(s)P (ex- 
cluding the p — UJ interference region) with five parame- 
ters: the magnitude and phase of the p' and p" contri- 
butions to Fk and Cj'^'^ -the magnitude of the qq con- 



tinuum, Eq. (11). In Fig. [sj we display our results for 
the electromagnetic pion form factor in the range - 
40 GeV^ < s < 10 GeV^. We confront them with the 
available experimental data for the electromagnetic form 
factor for -10 GeV^ < s < 10 GeV^ and the transition 
form factor for -40 GeV^ < s < -0.8 GeV^ (both are 
normalized to 1 at s = 0). First, we note that in the 
space-like region the data sets for the two form factors 
look identical (this is not expected to be the case for 
the time-like region since, for example the a;(782) only 
contributes to the transition form factor). One can see 
that our model describes all the available data through- 
out the shown kinematics. This serves as an a posteriori 
evidence that this s-dependence is in both cases driven 
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FIG. 3; Our results for the pion electromagnetic form factor 
for = 1 GeV^ (solid line) and ^? = 10 GeV^ (dashed line) 
vs. experimental data on the time-like and space-like e.-m. 
form factor from [25] (solid circles). 

by the same mechanism. In the case of the electromag- 
netic form factor, our result is a prediction for the s- 
dependence at large |s|, where no data exist so far. In 
particular, we predict that, as for the transition form 



factor case, \sF2t^(s)\ has to rise asymptotically roughly 
as s^/'^, unhke pQCD predictions that feature at most a 
logarithmic limit for that combination. 

To summarize, we presented a calculation of the tran- 
sition and electromagnetic form factors of the pion. We 
used dispersion relations to provide a unified description 
of the pion form factors in the time-like and space-like 
regions. In the hadronic energy range, we accounted 
for hadronic, resonance mechanisms in a fully unitarized 
manner. For asymptotic contributions, we proposed a 
new mechanism that features a reggeized quark exchange. 
We relate the parameters of such an exchange to tttt scat- 
tering data and deduce that the quark-Regge intercept is 
approximately a<j(0) '-^ 3/4. Using this value as input for 
the asymptotic behavior of the pion form factors, we ob- 
tain for both sF{s) oc s^/^, in agreement with the recent 
BaBar data. We notice that when the transition form 
factor is renormalized so that ^^-,,(0) = 1 its dependence 
on s in the space-like region is consistent with that of the 
electromagnetic form factor, as shown by the open circles 
in FigjSj We use the normalization of the Regge-behaved 
tqq,2TT and tqq^TTj as the only free parameter, and are able 
to describe all available data on pion form factors. 
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